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. ' Abstract 

C^ ' An energy-momentum tensor for general relativistic spinning fluids compatible with Tulczyjew-type 

supplementary condition is derived from the variation of a general Lagrangian with unspecified explicit 
form. This tensor is the sum of a term containing the Belinfante-Rosenfeld tensor and a modified perfect- 
OO ' fluid energy-momentum tensor in which the four-velocity is replaced by a unit four-vector in the direction 

of fluid momentum. The equations of motion are obtained and it is shown that they admit a Friedmann- 
Robertson- Walker space-time as a solution. 
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1 Introduction 

According to the general theory of relativity, free particles move along trajectories which are geodesies of the 
background space-time. Kinematically, the trajectory of each particle can be characterised by the particle 

f^ I four-velocity which is a time-like vector tangent to its path. This tangent also shows the direction of the 

particle momentum. On the other hand the dynamics of a general relativistic perfect fluid can be described by 

CO I its energy-momentum tensor. This energy-momentum tensor is expressed in terms of the fluid energy density 

l/^ ■ and pressure and a projection tensor which is itself constructed out of the four-velocity of the flow-lines and 

^^ I the space-time metric. 

^^ ■ When the particle has an internal classical spin, it no longer moves along a geodesic path in general due 

to a coupling between the particle spin and the space-time curvature. The dynamics of a spinning particle 
is described by the Mathisson-Papapetrou-Dixon (MPD) equations [T], according to which the particle four- 
momentum is time-like but not co-linear with its four-velocity. Thus in treating the dynamics of a spinning 
particle, we deal with two distinct time-like four-vectors corresponding to the velocity and momentum foul- 
ed [ vectors. These two preferred time-like vectors have been used to construct the two most favoured supplementary 
conditions (needed to render the MPD equations complete), the so-called Pirani condition and the Tulczyjew 
condition, see e.g. [T]. In the former the particle spin tensor has no electric components in the instantaneous 
zero-velocity frame while in the latter it has no electric components in the instantaneous zero-momentum 
frame. 

Several extensions of perfect fluid to spinning fluids, i.e. perfect fluids whose particles have an internal but 
classical spin, have been made in the literature basically by adding terms containing spin to the perfect-fluid 
energy-momentum tensor [2] El |4l [5] . The interrelations between the dynamics of spinning particles and the 
dynamics of spinning fluids have been studied in [B] in the framework of Einstein-Cartan theory. In these 
studies always a Pirani-type supplementary condition has been used. The question then arises if a formulation 
of spinning fluids consistent with Tulczyjew condition is possible. The aim of the present work is to address 
this question. We show that such a formulation is possible provided the projection tensor which constitutes 
the energy-momentum tensor, is constructed out of a time- like four- vector in the direction of the momentum, 
rather than the velocity four-vector. 
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To construct such a spinning fluid model, we deploy the machinery developed in ref. [S] in which a spinning 
fluid energy-momentum has been derived from a general Lagrangian whose explicit form is not specified. We 
start with such a Lagrangian and end up with an energy-momentum tensor whose (explicitly) spin-dependent 
terms are similar to the corresponding terms in the above mentioned reference but in the remaining part the 
usual perfect-fluid projection tensor is replaced with a modified one. In this modified projection tensor the 
four- velocity is replaced with a unit four-vector in the direction of the four-momentum. 

In the following sections we first briefly review the MPD equations describing the motion of spinning particles 
in curved space-times. Then we follow a customised version of the procedure developed in [5] to obtain the 
energy-momentum tensor. The equations of motion can then be obtained by imposing the requirements of 
symmetry and conservation on the energy-momentum tensor. We use the resulting equations of motion to 
study the dynamics of the spinning fluid in a flat Friedman- Robertson- Walker universe generated by the fluid. 
We show that these equations admit a solution of this type. A discussion of the results concludes the main 
sections. Some supplementary calculations are collected in an appendix. 

2 The MPD Equations 

The motion of a spinning particle moving in a curved space-time is described by the so-called MPD equations. 
Neglecting the particle multi-poles higher than dipole, these are given by [T] 

P^ = -^i^^a/S^-'s^^ (2) 

where s'^'^ are the particle spin tensor, p^ are the particle momentum, R^j^^nX represent the curvature tensor, 
over-dots mean covariant differentiation ^"Vq., and v^^ are the particle velocity. These equations supplemented 
by the so-called Tulczyjew condition 

P^s^'" - (3) 

result in 

PpP^ — const. — —m , 
Sij.uS^" = const. = 2s^ 

where to, s are the particles mass and spin respectively. Hence, defining 

m 
the normalization 

is guaranteed by the equations of motion. The 4- vector x'^, sometimes called the particle dynamical velocity, 
is not in general co-linear with the particle kinematical velocity, u'' . 

There is no equation of motion for w'^ and it should be determined indirectly from the other equations of 
motion. These equations result in the following relation [71 [8] 

which simplifies further when a specific gauge, say 

fKX" = -1 (5) 

is chosen. This is the gauge in which the instantaneous zero-momentum and the zero-velocity frames are 
simultaneous pj. The MPD equations do not guarantee that v^Vf^ is constant. 



3 The energy- momentum tensor 

We construct a spinning fluid model in which the fluid equations of motion reduce to the MPD equations for 
a single spinning particle. We begin with a matter Lagrangian of the form 

L = LK,e^',e^',X\a^X') (6) 

where 

m'' = — — 
dt 

are the fluid four-velocity, X^ are the identity coordinates, and e^ are the tetrad fields. The later satisfy 

Gai^e'^; = rjab, (7) 

a, b ranging from to 3. The fields X^ , i = 1, 2, 3 label the flow-lines in the sense that they identify the flow- 
line passing through any given point of space-time. To generate these, one can choose an arbitrary space-like 
hyper-surface equipped with a coordinate system X^ and label each flow-line by the coordinate value of the 
point where it intersects that hyper-surface [3]. 
We require 
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d^^X' = (8) 



which restricts the flow-lines to be directed along the integral curves of X* — const. The entropy per particle 
is assumed to be conserved along these curves. 
The conjugate momentum is defined by 

^' = S- W 

and the spin (density) tensor by [5] 

dL 
where 1°. = — — n- . This may be thought of as the density of the sum of individual particles spin s^'^ . For later 

/^ oe^ 
convenience we also define 

TT^^ , (11) 

representing a unit four- vector along the fiuid momentum. 
The action 

S = - yf^Rd'^x + I yf^Ld'^x (12) 

describes the matter and the space-time dynamics. Requiring this action to be stationary under variation of 
tetrad field, and making use of the Einstein equation 

G^v — —KTf^i, (13) 

we obtain 

By noting that the variation of L with respect to the tetrad field gets contribution from both explicit dependence 
on the field and implicit dependence on it via covariant derivatives, we have 



where 



here the first term in the parentheses comes from variation of ^/— g, the second two terms from variational 
derivative with respect to tetrad and the last term from implicit tetrad dependence via connection coefficients, 
as is shown in the appendix. Terms like u^S°"^ are particular cases of the Belinfante-Rosenfeld tensor [2]. 
On the other hand, by using the transformation properties of L, it can be shown that [lOj 



dL , dL ., 


dL 


— i-e H i-e - 


u' 


del "^ del'^ 


du^ 



d{d^x 

Assume that the quantity d^X^ depends on the particle density number n. Thus L depends on n implicitly. So 
the last term in the above equation may be simplified further. Consider an infinitesimal flux tube cfiX whose 
sections satisfy equation ([8]). Let jij denote the metric of this subspace and 7 its determinant. This metric 
can be obtained from the space-time metric by 

Y'=g^''d^x'd,x^. (14) 

The number of particles in this tube is assumed conserved and is given by iV = n^p^. We have 

dL dLdn dj dL 



d{d^X') dnd-fdid^X') dn 

and hence 

dT dT 

——-d.X^ = n—j,,d'^X^d,X\ 
d[dfj,X^) dn 

We also have 

jijdf^X^d^X' = g^^ + Ufj,u^ (15) 

note that v^ is normal to the sections of the subspace. Thus we obtain 

^^ .d^X^ = ^nK (16) 



didt'Xi) '^ dn 



where /i^^ = 5^1/ + Ufj_Ui, is a projection tensor. 
Putting the above results together, we get 



dL _ . . dL, 1. 



where we have added a vanishing term 

A^^e'^'^Se^ = 

to ensure that we will end up with a proper energy-momentum tensor. The explicit form of yl^^ depends on the 
choice of tetrad. There are two natural choices for the tetrad field. These correspond to taking the time-like 
component of the tetrad along the directions of fluid four-velocity or four-momentum. (These choices have 
been used in [11] for the case of spinning particles in a different context). We choose a tetrad for which Cq 
coincides with the momentum direction, i.e., we set 

e'o = TT^ (17) 

In this case A^i,e°''^Se^ = is guaranteed by orthonormality of the tetrad components. It turns out that this 
choice is consistent with the supplementary condition 

P^S'"' = 0. (18) 

The other choice is given by Cq — u^ which will result in an energy-momentum tensor as the one given in [S] 
together with a different supplementary condition. 



Now let us define f^ as the difference between tlie two velocities u'^ and tt'^, 

f'' = u''- n". (19) 

This results in 

Ur + 2/p7r^ = 0. 

If we now neglect ffj.f'^, an assumption based on the expectation that the effect of spin is not too strong, we 
get 



Utt^' = 0. 



We now use this to obtain 



Uf.u^e'^'Se'', = {tt^, + U){7T. + U)e''^Se: 



^^.^.e^'^Je^ 



We therefore obtain 



where in the right hand side we have interchanged fi and v, and have defined 

(or, equivalently, H^j^^ = e^Cij,). Here P^ is given by relation ^. Also, the fiuid pressure is defined by 

P = L-n^. (20) 

Thus we reach at the following energy-momentum tensor 

T,,. = H^,P + P^u, - ^V^iS^^u, + S°,u^ + S^^u'^) (21) 

This is a generalization of the energy-momentum tensor given in 10 for spinning dusts by taking the fiuid 
pressure into account, and modifies the one given in [3] which rely on the Frenkel condition. It is also different 
from the energy-momentum tensor given in [S] by the form of the projection tensor. 

The second term in Eq. (pij) is the same as those of j3] and |S]. The third term is the divergence of some 
Belinfante- Rosenfeld tensors. In the first term the usual perfect fluid term g^^ + u^Ui, has been replaced by 

4 Equations of motion 

By requiring this energy-momentum tensor to be symmetric we obtain the equation of motion of the spin 

Pf,u, - P,u^ = V„(m" V) (22) 

Demanding now the energy-momentum tensor to satisfy the conservation equation VT^,^ — 0, yields 

Hf^.V'P + PV^ii:^!:,) + ^{P^u,) = -^R^.c.iiu" S'^^ (23) 

This equation can also be derived by variation of world-lines. 

By projecting Eq. (|23p parallel and normal to tt^, we can reduce it to the following equations 

PW^TT^ + W^{Eu^) = ii?^.a^^^u''^"^ (24) 



H^.V^'P + Ptt"V,it,, + Sm^V^tt^ = -iff^„i?\a/ju"5"'5 (25) 
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respectively. Here E denotes the norm of P^. The first of these is a conservation equation and the second 
is an equation of motion for the fluid momentum. 

For a spinning dust, P = and p — nm, m being the particles mass. Defining 



n n 



and using V ^{nu'-'') — 0, equations (P^ and (P^ reduce to Eqs. ([T} and ^ respectively, the MPD equations. 

5 Applications 

In this section we apply the above formalism to a flat cosmological model. Consider a flat Friedman-Robertson- 
Walker space-time with the following line element 

ds^ = -dt^ + a'^{t)S,jdx'dx^. (26) 

where i, j run over the spacial dimensions, a{t) is the scale factor. For co-moving matter we have m^ = (1, 0, 0, 0). 
Now as we do not expect the four-momentum to be co-linear with the four- velocity in general, we start with 
P^ = {p{t), ep^) where we will neglect e^ and higher order terms. Note that this satisfies u^^P^ = —p, p being 
the energy density. The later is the fluid analogue of ([5]). Thus we have E — p and 

1./ 

P 
By inserting this into the energy- momentum tensor (|2ip we obtain 

r°° = p, (27) 

T" ^ 4 (28) 

T°' = i (e{2w + l)p' ~ dtS'^' ~ — 5°A (29) 

where w is the so-called equation of state parameter given hy w = —. Now the Einstein's equation gives 



and 



^Oi I o'^(''-' cOi 



e{2w + l)p'it) = dtS'" + 2^S'>\ (32) 

a[t) 



On the other hand, the spin equation of motion (|22|) yields 



^ + 5«5-- = 0, (33) 

at a 



, 4-S'* = -ep\ (34) 

dt a 

The first of these immediately results in 

S'^t) ^ a'^ {a{t))-^ (35) 



in which tr*^ are integration constants. The second equation when combined with equation (I32p yields 

S°' = ea"''(a(t))-" (36) 

^ fTO'a(i)(a(t))-("+i) (37) 



where 



w+l 

Aw + 3 



w + l 
From the translational equation of motion (|25p we obtain 

e{l + w)^ + e{4 + 5w)-p' = --S"' (38) 

at a a 

for which we have aheady a solution, relations ([55]) and (1571) . Inserting this solution back into the differential 
equation shows that it is valid for any w j^ —1. One can easily check that the above solutions of the equations 
of motion satisfy the conservation equation (|24p however they are consistent with the supplementary condition 
([T8l) only for 

a°' = (39) 

i.e. we end up with a four-momentum parallel to the four-velocity in this case. This is due the maximal 
symmetries of the space-time described by (P5|. From equation ([55]) we have S^ — iS'^^S'"' = ^aija^^a{t))^^ . 
Thus the fluid spin (the spin density integrated over a spacial volume) is space-like and goes like {a{t))~^. 

For both the above model and the usual perfect-fluid model the application of the Einstein's equation gives 
the same expression for the energy density and pressure in terms of the scale factor. Thus, in the present 
context the evolution of a{t) is the same as the standard model with a usual perfect fluid. This is in agreement 
with the result obtained in [4] . However it would be possible to obtain a different evolution for the scale factor 
by applying an averaging on spins scheme, as it has been done in [T^] for the case of Einstein-Cartan theory. 

6 Discussion 

We have constructed a consistent formulation of spinning fluids compatible with the Tulczyjew-type supple- 
mentary condition by starting with a rather general Lagrangian depending on a tetrad field that its gyration 
simulates the spin. There are two natural choices for the time-like component of this tetrad. These correspond 
to the velocity and the momentum directions and result in different supplementary equations. Our choice 
advocate the supplementary equation Pfj_S^^ = 0. The resulting energy-momentum tensor reduces to the usual 
perfect fluid energy-momentum tensor in the case where the spin is turned off and to the energy-momentum 
tensor introduced in Ref.[10] for spinning dusts, and the relevant equations of motion reduce to the MPD 
equations. The general form of this energy-momentum tensor and consequent equations of motion differ from 
those of Ref. [5] mainly due to different forms of the projection tensors. We have shown that the equations 
of motion admit a flat FRW space-time as a solution. For this solution the four-velocity and four-momentum 
are co-linear due to maximal symmetries of the solution. The dynamics of space-time in this model is similar 
to the one obtained by a usual perfect fluid model but applying an averaging process could result in different 
evolutions. It would be interesting to find solutions for which the the fluid velocity and momentum are not 
co-linear. Research is underway along this direction. 
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A Tetrad variation of the Lagrangian 

In this appendix we calculate the variation of L with respect to tetrad. We have 



But 







= 


(^^'-'-(-H))-- 


= 




= 





and 
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2 V de% del, ^^l 

2 " V de% de% de^ 

rIT BT rIT \ 

1 / BT BT 



where we have omitted a surface term. Thus 



dL bi. 3i .h,. 1. 
Beb "" Beb 



6eL = ( ^e"- + j—e"- + ^W^iu'^S"'' + u'^S'^" + u^Sf^n ] eaJe^. 



References 



[1] 

[2] 
[3] 
[4] 
[5] 
[6] 

[7] 

[8] 

[9] 

[10] 

[11] 
[121 



W.G. Dixon, Proc. R. Soc. London, scr. A 314 (1970) 499. 

J.R. Ray, L.L. Smalley, Phys. Rev. D 26 (1982) 2619. 

J.R. Ray, L.L. Smalley, Phys. Rev. D 27 (1983) 1383. 

Y.N. Obukhov, O.B. Piskareva, Class. Quantum Grav. 6 (1989) L15. 

L Bailey, Ann. Phys. 119 (1979) 76. 

W. Kopczynski, Phys. Rev. D 34 (1986) 352. 

K.P. Tod, F. de Felice, M. Calvani, Nuovo Cimento, 34 B (1976) 365. 

J. Ehlers, E. Rudolph, Gen. Rel. Grav. 8 (1977) 197. 

J.D. Brown, Class. Quantum Grav. 10 (1993) 1579. 

L Bailey, W. Israel, Commun. math. Phys. 42 (1975) 65. 

G.E. Tauber, Int. J. Theo. Phys. 27 (1988) 335. 

M. Gasperini, Phys. Rev. Lett. 56 (1986) 2873. 



